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ABSTRACT 

The spectral line polarization of the radiation emerging from a magnetized astrophysical plasma depends on the state of the atoms 
within the medium, whose determination requires considering the interactions between the atoms and the magnetic field, between the 
atoms and photons (radiative transitions), and between the atoms and other material particles (collisional transitions). In applications 
within the framework of the multiterm model atom (which accounts for quantum interference between magnetic sublevels pertaining 
either to the same y-level or to different y-levels within the same teiTn) collisional processes are generally neglected when solving the 
master equation for the atomic density matrix. This is partly due to the lack of experimental data and/or of approximate theoretical 
expressions for calculating the collisional transfer and relaxation rates (in particular the rates for interference between sublevels 
pertaining to different y-levels, and the depolarizing rates due to elastic collisions). In this paper we formally define and investigate 
the transfer and relaxation rates due to isotropic inelastic and superelastic collisions that enter the statistical equilibrium equations for 
the atomic density matrix of a multiterm atom. Under the hypothesis that the interaction between the collider and the atom can be 
described by a dipolar operator, we provide expressions that relate the collisional rates for interference between different 7-levels to 
the usual collisional rates for /-level populations, for which experimental data or approximate theoretical expressions are generally 
available. We show that the rates for populations and interference within the same /-level reduce to those previously obtained for 
the multilevel model atom (where quantum interference is assumed to be present only between magnetic sublevels pertaining to any 
given /-level). Finally, we apply the general equations to the case of a two-term atom with unpolarized lower term, illustrating the 
impact of inelastic and superelastic collisions on the scattering line polarization through radiative transfer calculations in a slab of 
stellar atmospheric plasma anisotropically illuminated by the photospheric radiation field. 
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1. Introduction 



. . The intensity and polarization of the spectral line radiation emerging from an astrophysical plasma depends on the population and 
J> atomic polarization (i.e., population imbalances and quantum interference between different magnetic sublevels) of the lower and 
upper line levels at each spatial point along the line of sight (LOS). Determining the population and atomic polarization of such 
rS levels requires considering the interactions between the atoms and photons (radiative transitions) and between the atoms and other 
material particles, such as electrons, atoms, and ions (collisional transitions). This problem can be very complex, especially when 
" " it comes to modeling the spectral line polarization produced by the joint action of anisotropic radiation pumping and the Hanle and 
Zeeman effects in multilevel systems. 

Within the framework of the density-matrix theory of spectral line polarization described in the monograph by Landi Degl'Innocenti 
& Landolfi (2004; hereafter LL04), it is possible to develop a consistent set of equations for multilevel systems, either by neglecting 
(multilevel model atom) or considering (multiterm model atom) quantum interference between pairs of magnetic sublevels pertain- 
ing to different /-levels (with / the level's total angular momentum value). The relevant equations are the radiative transfer equation 
for the Stokes parameters (where the coefficients of the emission vector and of the propagation matrix depend on the values of the 
atomic density matrix) and the master equation for the atomic density matrix (which includes both radiative and collisional rates). 

While for the multilevel model atom LL04 derived the expressions for both radiative and collisional rates (assuming isotropic 
collisions), for the multiterm model atom they only provide the expressions for the radiative rates. The aim of this paper is to 
formally define the collisional rates for a multiterm atom, and to find their relevant properties, focusing our attention only on 
isotropic inelastic and superelastic collisions. The treatment of elastic collisions in a multiterm atom is actually more complicated, 
and will not be treated here. Such collisions (e.g., with neutral hydrogen atoms) tend to equalize the populations of the sublevels 
pertaining to any given /-level and to destroy any quantum interference between pairs of them. A similar depolarizing role may be 



Article number, page 1 ofll3l 



caused by inelastic and s uperelastic collis ions between the 7-levels pertaining to any given term, especially when such 7-levels are 
very close in energy (see 'Bommier'"2009', for the hydrogen case). For the sake of simplicity, the latter type of collision will also be 
neglected, the investigation being limited to inelastic and superelastic collisions between different terms. 

In the main body of this paper, we provide suitable expressions for the transfer and relaxation rates caused by isotropic in- 
elastic and superelastic collisions, taking the possibility of atomic polarization in all the terms of the model atom into account. 
Particular attention is given to the collisional transfer and relaxation rates for interference between magnetic sublevels pertaining to 
different 7-levels, the physical ingredient that cannot be accounted for with a multilevel model atom. Since there are basically no 
experimental data for such rates, we provide approximate expressions here that relate such rates to the usual collisional rates that 
describe transitions between different 7-levels (for which experimental data or theoretical expressions are generally available). As 
a consistency proof of our derivation, we show that the transfer and relaxation rates for populations and interference between pairs 
of magnetic sublevels pertaining to the same /-level reduce to those derived by LL04 for the multilevel atom case. 

In the last section we present an illustrative application of the theoretical scheme developed here. We consider a two-term 
atom with unpolarized lower term, and we show the sensitivity to the collisional rates of the linear polarization of the radiation 
emerging from a slab of given optical depth, located at a given height above the "surface" of a solar-Uke star, and illuminated by its 
photospheric radiation field. 



2. Transfer rate due to inelastic collisions 

We consider a multiterm atom (see Sect. 7.5 and 7.6 of LL04) in the absence of magnetic fields, and we describe it by means of 
the density matrix elements ppisiJM, J'M'), with J the total angular momentum, M its projection along the quantization axis, L 
the orbital angular momentum, S the spin, and /3 the electronic configuration. This atomic model accounts for quantum interference 
(or coherence) between pairs of magnetic sublevels pertaining either to the same 7-level or to different 7-levels of the same term 
(7-state interference). We also work using the multipole moments of the density matrix (or spherical statistical tensors), defined by 
the equation 

^"p^(y,/) = ;^(-i)'-«V2FTT( ^ ^ \pisLsiJM,j'M'). (1) 

MM' \ 

Although collisional processes can be very eflicient in coupling 7-levels pertaining to the same term, in this investigation we only 
consider collisional processes coupling populations and coherence pertaining to different terms. 

In a given, although arbitrary, reference system, transfer processes due to inelastic collisions contribute to the time evolution of 
a particular density matrix element according to the equation 

j^P/}Ls(JM,J'M')= ^ C,{J3LSJMJ'M',l3tLcSJtMtJ[M'i)pi3,L,s{JcMc,J[M'^) , (2) 

PcUJcMiJ'iM'i 

where C/ is the inelastic collision transfer rate and where the quantum numbers (fi^LfS ) denote any term having energy lower than 
the term (j6L5)01 In a new reference system, obtained from the old one by the rotation R, recalling the transformation law (see 
Eq. (3.95)ofLL04) 

[pi^,s iJM, rM')l^^ = Yj ^NMiRT Km' (R) [PI^LS iJN, J'N')l^^ , (3) 

NN' 

with D-'i^j^iR) the rotation matrices, and its inverse 

[ppLs{JM,J'M')\^^ ^Y,^iN(R)Oi'N'iRr [PPLS^JN,J'N')\^^ , (4) 

NN' 



we have 

^ [ppLS (JM, J'M')l^^ = Z Z ^'^M(Ry ^N'M'iR) ^Im^R) ^i',M'(Ry cms JNJ'N',I3cUS JfNcJ',N',) 

PcLtJcMtJ'fM'i \nN'NiN'i j 

X [pp,L,s{JiMt,J'iM',)\^^^ . (5) 

The assumption of isotropic collisions implies that all the quantization directions are equivalent, so that Eqs. (13) and Q have to be 
identical. It follows that the collisional rates must satisfy the relation 

C,(l3LSJMJ'M',/3eL(SJcM(J',M'f) = ^ Oif^iRf 0';,,f^,{R)D'l^^{R)D'l^^,{Rr C,^^^ . (6) 

NN'NiN', 



' We assume that there is no overlapping in energy among the various terms of the model atom under consideration. 
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After coupling through Eq. (lA.lOl i the rotation matrices and £)-(„^,(/?), as well as the rotation matrices 2)^^^ (iJ)* and 

^iv'M' using the complex conjugate of Eq. ( lA.llb on the ensuing expression, Eq. (|6]l takes the form 

Ci{J3LS JMJ' M' ,l3tLtS JeMtJ[M'c) = ^ CiifiLS JNJ'N' ,/3{LeS J(NeJ(N'f) 

X i-lf''^^^'''^''(2K + 1)(2K' + 1){2K" + 1) 

KK'K" 

r J K \ ( r J K \ ( Ji f^' \i Ji J'e ^' 

X I A?' -N P j \ M' -M Q j \ Ne -N'^ P' j\ M( -M; Q' 

X I p p> P" j( g Q' Q" j®P"G"(^)- 

As the righthand side of Eq. ^ must be independent of the rotation R, the index K" can only take the value K" = 0. This implies 
K = K' , P ^ -P' , and Q ^ -Q' . Using Eq. iKA\ . we obtain 

CiipLS JMf M' ,PeL(S JeMeJ'cM't) = ^ C,(J3LS JNJ'N' J{NeJ(Nf) Y^(2K+l)(-lf-^^^'-^''-''-Q 

I r J K\l r J K\l J, 4 ^ ]( ] (8) 

\ N' -N P j \ M' -M Q j \ N( -N'l -P j\ M( -M'^ -Q j ' 
2.1. Multipole components of the inelastic collision transfer rate 

Introducing the multipole components of the inelastic collision transfer rate, defined by the equatiorQ 



Cf\l3LSJJ',l3tLcSJ[J[) 



'/ + /' + ! 



X (-l)"""'"""'''! j^' J^v ^p][n' -Ne ^p]ci(PLSJNJ'N',l5eUSJeNeJ'tN't), (9) 

NN'N,N[ ^ I \ ( I 

and making use of Eqs. dA.ll l and ( IA.2I ). Eq. ([8]) can be written in the form 

2(2^+l)(i', X lYriPLSJJ',p,L,SJ,r,). (10) 



X 

K 



Substituting Eq. ( fTOl i into Eq. (|2]), and recalling the definition of the multipole moments of the density matrix (see Eq. ([T])), with the 
help of Eq. (IA.3I ). we find the following equation for the spherical statistical tensors 



^^"p§(7,/)= ^I '^jl'^Jll Cf\l3LSJJ\PcLiSJeJ'i)^''"'pl{Ji,J'd ■ (11) 

Taking the complex conjugate of Eq. (|2| and recalling that ppisiJM, J'M')* - ppisU'M' , JM), we have 

CiifiLS JMJ'M' , PiLcS JcMiJ'iM'iY = CiifiLS J'M' JM, PtLeS J'^M'^JtMe) , (12) 
and therefore, using Eqs. (lA.H and (IA.2| |. 

Cf'>(J3LSJJ',P[LtSJi;J'i)* = i^-iy+J'-J'-^'i cf^(J3LSJ'J,PcLcSJ'eJd ■ (13) 
Setting K -Q 'm Eq. (|9]l, and using Eq. (IA.4I I. we obtain 

Cf\/3LSJJ',/3eLcSJcJ'e) = 6jr 6j,r V Ci{j3LS JNJN,/3cLeS JcNcJgNc) , (14) 

where the transfer rate CjifiLS JNJN,/3(L(S J(NfJ(N() is the usual (inelastic) collisional rate for the transition from the lower 
magnetic sublevel \/3(L(S JfNf) to the upper magnetic sublevel \/3LS JN), generally indicated in the literature with the notation 
Ci{J3(L(S JcN[ — > pLS JN). Since this rate is non-negative, the 0-rank multipole component is also non-negative. 



^ The factor yJ + J^+l/ + 7^ + 1 is introduced in order to get simpler relations between these rates and the usual collisional rates connecting 
atomic populations. This factor reduces to the one introduced in the multilevel atom case (see Eq. (7.87) of LL04) when interference between 
different /-levels is neglected (/ = J' and Jc = J'(). 



Article number, page 3 ofll3l 



2.2. Relations with tine collisional rates for J -level populations 

In most cases, the only collisional rates for which experimental data, or approximate analytical expressions, are available are the 
collisional rates connecting the populations of different /-levels (following the notation generally used in the literature, these rates 
will be indicated through the symbols CiificLcS J( — » PuLuSJu) and CsifiuLuS Ju — » fi(L(SJi), the indices / and S standing for 
"inelastic" and "superelastic", respectively). It is important therefore to find suitable relations between such rates and the collisional 
rates introduced in this paper for a multiterm atom. 
Observing that 

CmitSJt ^ PuUSJu) = y CiiPdcSJcNc ^ PuLuSJuNu) , (15) 

from Eq. (fT4l i we immediately have 

cf(J3LSJJ,l3[LeSJcJe) = Ci(J3[LeSJc jiLSJ) . (16) 

In Eq. if we couple through Eq. ( lA.lOb the rotation matrices Dj^j^iR)* and Dj^^j^^^(R), as well as the rotation matrices £)^,^,(7?) 
and l/^,^,{Ry, by requiring that the ensuing expression is independent of the rotation /?, we find the relation 

CiipLS JMf M' ,l3fLiS JeMcJ'fM'i) = ^ C,(J3LS JNJ'N' ,(3eL{S JiNtJ'iNf) ^(2/: + 1) (-l)'^-*'+^^*'^^-2 



J Jc K \ l r J'l, K\l J J[ K \ l r J'f K 
-M M[ Q ]\ -M' M'^ Q ]\ -N N( P ]\ -N' N'^ P 



(17) 



Defining a different set of multipole components of the inelastic collision transfer rate through the equation 

Tf^{pLSJJ',lieL(SJ(J',)^j^^^^ 2 {-if'-'^'i^^ ?)(-"^' N' p ]ci(J3LS JNJ'N' ,/3eLcS JcNcJ'eN',) , 

(18) 

Eq. ( [TtT i can be written in the form 

Ci{pLSJMJ'M',l3eLeSJeMeJ'(M'f) ^{-\)^'r^' {J e + J't + ^)Yu[-M M( g ) ( -M' M', Q ^^(^1^^ J f ,l3cLcS JeJ'i) . 

(19) 

As pointed out in LL04 for the multilevel atom case, this decomposition of the collisional rate has an interesting physical interpreta- 
tion, because it shows that the interaction between the atomic system and the collider can be described by a sum of tensor operators 
of rank K acting on the state vectors of the atom. Starting from Eq. (|9) and using Eq. ( IA.8b . after some algebra the following relation 
between the multipole components C^^^ and V''^ can be found; 

Cf\pLSJJ',l3eLtSJcJ'e) = ^{J + J' + \)(Jc + J'f+ l)^(-l)^'+^^-'''+^ | ^,^Y^P(JiLS JJ' ,l3iLeS Jeff) . (20) 

For the K - Q multipole component, using Eq. (IA.7l l. we have (cf. Appendix A4 of LL04) 

cf^(J3LSJJ,f3cLeSJ(Je) = ^ vf^iJiLS JipdeS JfJd . (21) 

K 

When the interaction can be described through just one operator of rank K, then 

Cf^(fiLSJJ',l3tLtSJcJ[)^ ^(J + r + \){Jc + J'f + l)(-lY^'''-^^''^ j, ^^vf^ifiLSJJ^PcLtSJtJ'i) . (22) 

The multipole component of rank K of the diagonal rates ( J - J' and J( - J'f) is thus related to the multipole component of rank 
by the equation (cf. Appendix A4 of LL04) 

J J K \ i J J K 



Cf\liLSJJ,l3tLiSJeJd = ^ ' cfifiLS JiPtLtS JtJ() = {-lf \ ^ [ dliiUSJt ^ PLS J) . (23) 
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A similar relation for the nondiagonal rates (describing the transfer of 7-state interference due to inelastic collisions) cannot 
be obtained through symmetry arguments alone. Such a relation can, however, be derived if some simplifying hypotheses on the 
interaction between the atoms and perturbers are introduced. It is well known that in the case of electrons with much higher energy 
than the threshold energy (i.e. under the so-called Born approximation), the Hamiltonian describing the electron-atom interaction 
depends on the dynamical variables of the atom only through the dipole operator (a tensor of rank K - I). A collisional process in 
an optically allowed transition can thus be treated, in a first approxi mation, as a radiative transition, and the collisional rate can be 
expressed through the oscillator strength of the same transition (e. g. ISeatonll 1 9621 : 1 Van RegemorteJ 1 9621) . 

For more insight on the nondiagonal rates C^f^^ijSLS JJ' ,/3cL[S JeJ^), we assume that the electron-atom interaction is described 
by a dipolar operator, and we proceed by analogy with the multiterm atom radiative transfer rate due to absorption processes (T^). 
Setting K,- - (i.e. assuming an isotropic radiation field) in Eq. (7.45a) of LL04, we have 



TA(/3LSKQJJ',PcLeSKQJeJ',) =(2Lf + ^(27 + l)(2J' + l)(2Je + 1)(2J'^ + 1) 

BificLcS ^ I3LS)jI , (24) 





Jt 






Lt 










{ J'l 


r 






J 




{ J'^ 


r 


I] 



where 7° is the angle-averaged incident radiation field, and B(fi(LeS — > pLS) is the Einstein coeflicient for absorption from the 
lower term (J3(L(S) to the upper term (fiLS). We recall that this quantity is connected to the Einstein coefficients for the individual 
transitions between fine structure 7-levels of the multiplet by the relation (see Eq. (7.57a) of LL04) 

B(J3cLcS Jt ^ PLS J) = {2Le + l)(2y + D | ^ | BifitUS ^ /3LS) , (25) 



which implies (using Eq. ([A^ 

B(J3eLeS ^ /?L5 ) = ^ B(J3cLcS Jt ^ jSLS J) . (26) 
J 

By analogy with Eq. ( l26l l, we define an inelastic collisional rate for the transition from the lower to the upper term Ci(J3fLeS 
PLS ) through the equation 

CifieLiS ^I3LS )^Yj CiQScLeS Je ^ PLS J) , (27) 

J 

where the sum is extended to all the 7-levels of the upper term to which a given /-level of the lower term can be connected through 
an electric dipole transition. By analogy with Eq. (l24t . and taking the multiplying factor introduced in Eq. (|9]l into account (see 
footnote 2), we can write 



Cf\j3LSJJ',/3eLfSJfJf) =(2Lf + l)(-l)i+^'+-';+'^ , 



(J + J' + 1)(2J + 1)(2J' + l)(2Je + 1)(2/; + 1) 



Je + J' + l 



X 



r lc}{j, 7 l}{ 7; r I ^/sls) . (28) 



This equation can be used to calculate the multipole components of the inelastic collision transfer rates for 7-state interference from 
the values of the usual inelastic collisional rates for 7-level populations. As a proof of the consistency of Eq. (l28T l, we observe that 
the 0-rank multipole component is given by 

Cf\pLSJJ,li[LfSJtJc) = CiiPcLfSJt ^ PLSJ) = {2Lc + 1)(27 + 1) | M CAficLtS ^ pLS) , (29) 



which is the analogous to Eq. (|251 ). while using Eqs. ( IA.7l i and ( |29] |, the diagonal terms are given by 

J J K 

dP(J3LSJJ,licLtSJf Ji) = (-1)^ \ ^ CiiPcLfSJe /3LSJ) , (30) 

Je Jc 1 

which corresponds to Eq. (l23T l for K = I. 



The sum appearing on the righthand side of Eq. i26) does not depend on the particular /-level of the lower term that is considered. 
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3. Transfer rate due to superelastic collisions 

A similar reasoning can be followed for the transfer rates due to superelastic collisions. These transfer processes contribute to the 
time evolution of a particular density matrix element according to the equation 

j^PpLs{JM,J'M')^ 2 Cs(/3LSJMJ'M\p,L,SJuM,J',M'Jpp^,L,,s(J.M,,JlM'J , (31) 

I3„l„j„m„j;,m;, 

where Cs is the superelastic collision transfer rate and where the quantum numbers (fiuLuS) denote any term having energy higher 
than the term (J3LS). Following the same steps as in Sect.|2] it can be shown that under the assumption of isotropic collisions, the 
transfer rate Cs can be written in the form 



CsiPLS JMf M' ,PuLuS JuMuKM'^) = ^/ (-1) 



Z(2^+1)(m' -M -k l)cf\mJJ',l5.,USJJ'X (32) 



where the multipole components of the superelastic collision transfer rate, C^P , are defined by the equation 



Cf^{J3LSJJ',l3uLuSJJ',) 



J + J' + l 



X 2 iJ' -N p](n' -k, p]cs(J3LSJNJ'N',fi.,L„SJ„N,J:.K)- (33) 

NN'N„N'„ * I \ u " I 

By substituting Eq. ( l32l i into Eq. ( |3T] ). and recalling Eq. ([1), we find the following equation for the spherical statistical tensors 



A/^"p^(7,7')= X ^ TJj'Vl Cf(fiLSJJ'A,USJj:y'''''YQ(JuJ:) ■ (34) 
The 0-rank multipole component is given by 

Cf(J3LSJJ',PuUSJuJ'u) = Sir 6j,x E CsijSLS JNJN,/3,L,S JuNJM = 6jr Sjj^ CsifiuUSJu ^ PLS J) , (35) 

^ NN„ 

where Cs (PuLuS Ju — > PLS J) is the usual superelastic coUisional rate for the transition from the upper level | PuLuS Ju) to the lower 
level \/5LS J). When the interaction between the atomic system and the colliders can be described by means of a single operator 
of rank K, it can be shown that the multipole components of rank K of the diagonal rates {J = J' and J„ - J^) are related to the 
multipole components of rank by the equation 

J J K ] ( J J K 



Cf\pLSJJ,l3,L,SJuJu) = (-1)'' !:'^" ^" ^ ' cfipLSJjPuUSJuJu) = {-lf \ ^" ^" ^ { CsiPuLuSJu ^/3LSJ) . (36) 



Ju Ju ^ j \ Ju Ju ^ 

As discussed in the previous section for the case of inelastic collisions, a similar relation for the nondiagonal rates (describing the 
transfer of 7-state interference due to superelastic collisions) can be derived under the assumption that the electron-atom interaction 
is described by a dipolar operator By analogy with the expression of the multiterm atom radiative transfer rate due to stimulated 
emission processes (Tj , see Eq. (7.45c) of LL04) in the presence of an isotropic incident field, we find the following relation 



Cf^(pLSJJ',j3„LuSJuJ'u) =(2L„ + 1)(-1)' 



(J + J' + 1){2J + l)(2J' + l)(2Ju + 1){2J;, + 1) 
Ju + J'u + l 



J Ju 1 
J'„ J' K 



Ju 7 l]{ru r I ]csiJ3uLuS^/3LS), (37) 

where we have introduced the superelastic collisional rate for the transition from the upper to the lower term Cs (fiuLuS — > /3LS ), 
defined by 

Cs (JiuLuS -^pLS) = Y,^s ifiuLuS Ju ^ PLS J) , (38) 
J 

the sum being extended to all the /-levels of the lower term to which a given 7-level of the upper term can be connected through an 
electric dipole transition^ 



The sum appearing on the righthand side of Eq. ( I38t does not depend on the particular /-level of the upper term that is considered. 
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4. Relaxation rates due to inelastic and superelastic collisions 

In a given reference system, relaxation processes due to inelastic and superelastic collisions contribute to the time evolution of a 
particular density-matrix element via an equation of the form 



jPpLs{JM,J'M') = - ^ [f(J5LSJMJ'M'J''M'')ppis{JMJ''M'') + g(J5LSJMJ'M'J''M'')ppis{J''M'',J'M')\ 

J"M" 

The conjugation property of the density-matrix elements (p/jts {JM, J'M')* - ppis (J'M', JM)) requires that 

g(J3LSJMJ'M'J"M") = f(J3LSJ'M'JMJ"M"y , 

so that Eq. i39[ can be written in the form 

d 
df 



PpLs{JM,J'M')^- ^ 



^5 (JiLS JMJ'M'f'M") ppLs (JM, J"M") + ^5 (JiLS J'M'JMf'M")* ppis {J"M", J'M') 



(39) 



(40) 



(41) 



dt 



In a new reference system, obtained from the old one by the rotation R, recalling Eqs. Q and (|4]i, we have 
[pi,Ls{JM,J'M')]^^^ = - 2 1^ Z ^NMiRy^N'M'iR)S(13LSJNJ'N'J"N") 

J"M"M"' ^ NN'N" 

+ \ Tj ^iM(RyK'M'(R)S(/^LSJ'N'JNJ"N"y 

NN'N" 

X Di;:,^„(R)K^,„iRy [ppls(j"m",j'm"')1^ 

Due to the isotropy of collisions, Eqs. (HTb and ( l42l i must be identical, which implies 
SipLS JMJ'M'f'M") 6 MM'" = Yj ^NM(Ry^N'M'(R)S(JiLSJNJ'N'J''N'')DiM,,,(R)D'!!,:,j^,,{Ry, 



(42) 



(43) 



NN'N" 



regardless of the rotation R. This requires the rate S (J3LS JNJ' N' J" N") to be independent of the quantum number (if not, the 
righthand side of Eq. (l43T l would not be zero for M + M'" , no matter the rotation /?). We can thus carry out the summation over 
via Eq. ( |A.9t to get (with the help of Eq. ( lA.lOl )) 





J" 




{ '\ 


J" 




['^' 


-N" 




\ M' 


-M" 


I) 



S(J5LSJJ'M'J"M") = Yj S(/3LSJJ'N'J"N") (-1)'^"-^" ^P-K + 1) ( ^' J 

N'N" K ^ 

Since the righthand side of Eq. (l44l l must be independent of the rotation R, index K can only take the value K - 0, which implies 
g = P = 0, A^' = A^", M' = M", and /' = J" from Eq. (IaH i. We thus obtain 

S{J5LSJJ'M'J"M") = 5m'M"6j'J" ^ ^ ^ ^ (y3L^ //A^ V'A^O ■ (45) 



Substitution into Eq. ( |4TI) gives 
d 



df 



ppLsiJMJ'M') = -So(/3LSJJ')pfiLs(JM,J'M') 



or, in the spherical statistical tensor representation, 

f/'VoiJJ') = -So{/3LSJJ')'''Yq(JJ') , 

where we have introduced the collisional relaxation rate 



So(/3LSJJ') = - 



YS (JiLS J J'M' J'M') + Yj^(I3LSJ'JMJM)* 

M' M 

The diagonal element 

Soi^LSJJ) = JJ—^^Y [S(/3LSJJMJM) + S {jSLS J J M J M)*] = ^j^Re ^ 5 08L5 77M7M)) 



(46) 
(47) 
(48) 

(49) 
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coincides with the one defined in LL04 for the case of a muhilevel atom. 

As shown in LL04, the diagonal element S o(J3LS JJ), which represents the relaxation rate of populations and interference 
between magnetic sublevels pertaining to the same 7-level (see Eqs. (l46l l and (l47b ). is connected to the 0-rank multipole components 
of the inelastic and superelastic collision transfer rates by the equation 

So(J3LSJJ)^ 2 cf^(J3,L,SJJ„/3LSJJ)+ ^ cf^(J3eLeS JeJe,/3LS JJ) . (50) 

To obtain a similar relation for the nondiagonal elements {S oifiLS J J') with J + J'), which represent the relaxation rate of J- 
state interference due to inelastic and superelastic collisions, we make the assumption, also in this case, that the interaction between 
the atoms and colliders is described by a dipolar operator, and we proceed by analogy with the multiterm atom radiative relaxation 
rates due to absorption (R^) and stimulated emission (Rj) processes (see Eqs. (7.46a) and (7.46c) of LL04). Assuming an isotropic 
incident radiation field (i.e. setting K, = 0), such radiative rates assume the simple form 

Ra(J3LS KQJJ'KQJf ) = ^ BiJiLS ^ l3„h,S )JI , (51) 

P„L„ 

Rs(/3LS KQJJ'KQJf) = ^ B(/3LS /3(LeS)J^ . (52) 

l3rL[ 

Introducing the inelastic and superelastic collisional rates for transitions between different terms (see Eqs. (l27T l and (l3Ft ). we have 
So(J3LSJJ') = ^''^P^^ ^ A<i«5) + X Cs(j6i5 ^ PeLiS) 

= Ci(J3LS J ^ /3,L,S 7„) + Y'^s(J3LSJ^ /SfLfS J() = S o(J3LS J J) = So(J3LS). (53) 

l3„L,iJ„ PrLiJr 

The relaxation rate of 7-state interference due to inelastic and superelastic collisions thus coincides with the relaxation rate of J- 
level populations and of interference between magnetic sublevels pertaining to the same 7-level. This rate, on the other hand, does 
not depend on the quantum number J, and is thus identical for all the /-levels of a given term. 

When collected together transfer and relaxation rates, the statistical equilibrium equations for the spherical statistical tensors 
can be written in the form 

^ ^"p^(7, J') = Y VttJtT '^f '^P^^ JJ'^PcLtS JtJ'e) ^'""'pliJt, J'l) 

PfLfJfJ'^ 

+ Z V '^jIj'Vi JJ'^/3uL,S JJ:,) ^"^"V§(7„, J'J 

-So(J5LSJjy^%(JJ') . (54) 



5. Application to the case of a two-term atom with unpolarized lower term 

We consider a two-term atom and denote the quantum numbers characterizing the lower and upper term by i/3(LfS ) and (J3uLuS ), 
respectively. The time evolution of the spherical statistical tensors of the upper term, when taking both radiative (see Eq. (10.1 15) 
of LL04) and collisional (inelastic and superelastic collisions only) processes into account is described by the equation 

^P'^'Yq(J„ j;j ^-2m Y Np„LAKQJuJ[r K'Q'J'Jj;,") ''"'"■Vq,(J:! , J',!') 

K'Q'JL'JL" 



+ Y '^A{liuLuSKQJ,J[„l3tLiSK'Q'JtJ'f)P''^''pl,{Je,J'c) 
Y ^e(J3uLuSKQJJ[^K'Q'j:j:') 



K'Q'J"J'„' 



+ Rs(/3,h,SKQJXK'Q'J','J::') 



B„L„S^K',jii Jill 



- So(fiuL„SJj:,) ^""-"'pliJu, J',) , (55) 

where is the magnetic kernel (see Eq. (7.41) of LL04), the radiative transfer rate due to absorption, while R^ and Rj are the 
radiative relaxation rates due to spontaneous and stimulated emission, respectively. 
We now make the following simplifying assumptions: 
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- There is no magnetic field. Under this assumption the kernel takes the simpler form 

Np„L.,s {KQJJl K' Q'J'JO = 6kk' Sqq. 5j^,r; brj- vp,,L..sj.„p.,L„sr., , (56) 

with Vjj^^L^^s j,„j3„L„s 3[, - [E(J3uLuS Ju) - E(J3uLuS J'J]/h, where E(J3LS J) is the energy of a given fine-structure 7-level, and h is the 
Planck constant. 

- The radiation field is weak so that stimulated emission can be neglected (Rs - 0). 

- The lower term is unpolarized (i.e., the magnetic sublevels of the lower term are evenly populated and no interference is present 
between them). Under this assumption the spherical statistical tensors of the lower term are given by 



Pe(y,,7,) = 5^o%,5.,.;^25 + l)(2L,+ l)77' ^^'^ 

where N is total number density of atoms, and the number density of atoms in the lower term. 
- The electron-atom interaction is described by a dipolar operator Under this assumption, defining through Eq. (l3Ft a superelastic 
collisional rate for the transition from the upper to the lower term (Cs (fiuLuS — » /3(L(S )), the collisional relaxation rate is given 
by (see Eq. Q) 

So(/3uh,SJX) = CsifiuLuS ^ l3eLrS) . (58) 

Taking the above-mentioned assumptions into account, andrecalhng that (seeEq. (7.46b) of LL04) ReQ^uLuS KQJuJ'uK' Q' J'^^ J'J') - 
6kk' Sqq' 6j,x 6j:x"A(J3yLuS /3(L(S), we obtain 

^/'''"VQiJu, = - 27riv;s„L„sy,„A,L„sy,', ''"'-"Yq(Ju, J'J 

ZV27f + 1 Nf 
T,(/S.L^SKQJ..JMSOOJM (2,,l)(2L,.l) 77 

- A(ji,us ^ ptLiS) pI{j„ 

V27f + 1 Nf 



\){2Lc+l)N 

- Cs(J3uL,S ^ j3(LeS ) ^"^"^p^(7,„ J'J . (59) 

As expected, under the hypotheses of isotropic collisions and unpolarized lower term, transfer processes due to inelastic collisions 
only contribute to the time evolution of the 0-rank spherical statistical tensors of the upper term. Assuming that the colliding particles 
are characterized by a Maxwelhan velocity distribution, the collisional rates Cf^ and cf^ can be related through the Milne-Einstein 
relation 

'E(/3„L„SJ,)-E(j3eLcSJe) 



KbT 

where T is the electron temperature. Using Eq. (l3Ft . the fourth term on the righthand side of Eq. ( |59] l can be written in the form 



(60) 



where Bt(vo) is the Planck function in the Wien limit (consistently with the assumption of neglecting stimulated emission) at 
temperature T, and where vq is the Bohr frequency corresponding to the energy difference between the centers of gravity of the two 
terms. (We neglect the frequency differences among the various components of the multiplet in the exponential appearing in the 
Milne-Einstein relation.) 

Taking the expression of T a(J3uL,iS K QJu J ^i,/3eLeS 00 J(Je) (see Eq. (10.124) of LL04) into account and performing the sum over 
J( using Eq. (IA.6I ), the second term on the righthand side of Eq. (|59l ) is given by 

We recall that the quantum theory of polarization described in LL04 is valid under the so-called flat spectrum approximation (that 
is, the incident radiation field that produces optical pumping in the atomic system must be flat over a frequency interval Av larger 
than the natural width of the atomic levels, and, when coherence between nondegenerate levels is involved, Av must then be larger 
than the corresponding Bohr frequency). For this reason, it is sufficient to evaluate the radiation field tensor Jq (see Eq. (5.157) of 
LL04 for its definition) at a single frequency within the multiplet. 

In stationary situations, recalling the relations among the Einstein coefficients 

BiPeLeS ^ /J„L„5) = ^^^B(/3„L„5 /3(LeS) = l^" ^ ] J", , AiJi^US ^ pdeS) , (63) 
IL( + 1 IL( + 1 2«Vq 
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the spherical statistical tensors of the upper term are given by 



Pq( u, u) + 1)(2L„ + 1) ^ A(/3„L„5 ^ PcUS) 



^ 1 + e' +27riy^„i„sy„,;3„/.„sy;,/A0S„L„5 ^ PiLtS) 

where, in analogy with the two-level atom case, we have introduced the quantity 



(64) 



e . (65) 

A{p,U,S PcLcS) 



It can be easily proven that if 5 =0, so that the upper and lower terms are composed by a single fine-structure 7-level, the expression 
of a two-level atom is recovered (see Eq. (10.50) of LL04, with //„ = = 0). 

Substituting Eq. ( l64l l into Eq. (10.127) of LL04, and introducing the frequency-integrated absorption coefficient of the multiplet 

^ = ^MeBiPeLeS ^ ^„L„5 ) , (66) 

and the absorption profile of the multiplet (in the absence of magnetic fields, and for the case of a two-term atom with unpolarized 
lower term) 

(2Je + l)(2J, + 1) { Lu L( 1 



f^^'^-L \ Jf Ju 5 ^^^AA^^-'/'.^.^y.-v), (67) 



where (^(vo - v) are Lorentzian profiles centered at the frequencies of the various components of the multiplet, we find the following 
expression of the emission coefficient in the four Stokes parameters: 

e,(v,n) X Z (-l)^-^^^^"^^"^^'^'^^^3(27„ + 1)(27: + 1)(27, + 1) 

KQ J«J'„Jc 

U Li I \ { U L( \ \ { I \ K\{ \ I K \ { U U K 
J( 7„ S \\ Jc J'u S \\ 7„ J:, Jt \\ K L„ Lf \\ J, S 



n-Kc n\ jK , ^ 1 ^(Vfi„L„S J„,p,L[S J, - v) + ^(Vfi^,L„S j;,.li,LrS J, - v)* 
X J riKl,^!-) J n(yn) 



■l4,BT(vo)<p(v)6ifi, (68) 
with / = 0, 1, 2, and 3, standing for Stokes /, Q, U and V, respectively. Here, v and Q are the frequency and propag ation direction 



l+e' 
, 2, anc 

of the emitted radiation, respectively, '7~^(/, Q) is the geometrical tensor introduced bv lLandi Degrinnocentiin983h . and <l)(vo - v) 



I propag a 
l ([T983h . 

is the complex emission profile 



(D(vo - v) = (p(vo -v) + i if/(vo - v) , (69) 

with 0(vo - v) a Lorentzian profile and i/'( vq - v) the associated dispersion profile^ The last term on the righthand side of Eq. (l68t rep- 
resents the contribution to the emission coefficient coming from collisionally excited atoms. Since collisions are assumed isotropic, 
this term only contributes to Stokes-/. 

As an example suitable to illustrating the sensitivity of the emergent scattering line polarization to the studied collisional rates, 
we consider a constant-property slab of stellar atmospheric plasma located at a given height above the surface of a solar-like star 
and characterized by a given optical depth Ar. Neglecting limb-darkening effects, the radiation illuminating the slab from below is 
characterized by an anisotropy factor w - ^/iJ^/j''^ given by 

cosQ'(l -I- cose) 

w = — , (70) 

where a is half the angle subtended by the stellar disk, as seen from the slab. We solve the equations of the non-LTE problem 
described in this section for the case of a ^S - transition, using the level energies and transition probabilities of the Mg ii li and 
k lines, and we calculate the fractional linear polarization of the radiation emerging at fi - cos 6 = 0.1, with 9 the angle formed by 
the local vertical (perpendicular to t he slab) and th e emission direction. The non-LTE radiative transfer problem is solved following 
the numerical methods described in iTVuiili o Bueno & Manso Sainz (1999). 



^ The equations derived here are valid in the atom rest frame. Nevertheless, under the assumption of complete redistribution on velocities (see 
Chapter 13 of LL04), the same equations can also be applied in the observer's frame, with 0(vo - v) and ifi(v(, - v) the Voigt profile and the 
Faraday- Voigt profile, respectively (provided that the atoms have a Maxwellian distribution of velocities). 
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Fig. 1. Left panel: QII profile of the radiation emitted across a - "P transition (the level energies and transition probabilities are those of 
the Mg II h and k lines) as obtained for dilferent values of the parameter e' (indicated in the plot). The arrows point to the wavelength positions 
of the two lines. Right panel: zoom of the line-core region of the 1/2 - 3/2 transition. We consider the radiation emitted at yu = 0. 1 by a slab 
located 0.03 stellar radii above the surface, and with an optical depth (at the line-center frequency of the 1 /2 - 1/2 transition) At = 0.5. We solve 
the full non-LTE radiative transfer problem within the slab, the boundary condition being the stellar radiation illuminating the slab from below 
(limb-darkening elfects are neglected). We consider a Doppler width of 26 mA, corresponding to a temperature of lO'' K and a microturbulent 
velocity of 1 km/s. We include the effect of an unpolarized continuum characterized by an opacity 10* times less than the line opacity at the 
line-center frequency of the 1/2 - 1/2 transition. The reference direction for positive Q is the parallel to the closest limb. 



Figure 1 shows the fractional linear polarization Q/^ pattern calculated for different values of e', considering a slab located 0.03 
stellar radii above the surface (corresponding to about 2 x 10^ km in the solar case), and characterized by an optical depth (at the 
line-center frequency of the 1 /2 - 1/2 transition) At = 0.5. We assume a Doppler width of 26 mA, corresponding to a temperature 
of 10^ K, and a microturbulent velocity of 1 km/s. The damping constant is consistently calculated as 



Avd A-kAvd 

withy,, = AifiuLuS /3eLfS)+Cs(fiiiL,iS jSfLf 5) the inverse life time of the upper term (elastic collisions, hence their broadening 
effect, are neglected). We include the contribution of an unpolarized continuum characterized by an opacity rj'j = 10"^ '//(vi/2-1/2), 
with Tj^jivi /2-\/2) = ^p(vi/2-i/2) the line opacity at the line-center frequency of the 1 /2 - 1/2 transition. We first note that in the 
slab model that we have considered (in which radiative transfer effects are significant), the Q// profiles show the typical signatures 
of y-state interference, such as the sign-reversal between the two lines, and the high polarization values in the far wings (see Stenflo 
1980, LL04, and Belluzzi & Trujillo Bueno 201 1). As can be observed, the modification of the scattering line polarization pattern 
due to inelastic and superelastic collisions (quenching effect) becomes appreciable only for rather high values of the collisional rates 
(e'>10-2). 



6. Conclusions 

In this paper we have formally defined and investigated the transfer and relaxation rates due to isotropic inelastic and superelastic 
collisions that enter the statistical equilibrium equations for the atomic density matrix of a multiterm atom (i.e., a model atom 
accounting for quantum interference between magnetic sublevels pertaining either to the same 7-level, or to different /-levels 
within the same term). 

While the numerical values of the collisional rates for 7-level populations are generally available (either from approximate 
theoretical expressions or form experimental data), the values of the collisional rates describing the transfer and relaxation of 
quantum coherence are in most cases unknown. In this work we focused our attention on the collisional rates for 7-state interference 
(the physical aspect that cannot be accounted for with a multilevel model atom). Under the assumption that the interaction between 
the atom and the perturber is described by a dipolar operator, we derived suitable relations between such rates and the usual 
collisional rates for 7-level populations. In particular, we showed that the collisional relaxation rate for /-state interference coincide 
with the relaxation rate for 7-level populations and for interference between magnetic sublevels pertaining to the same 7-level. We 
also observed that this rate does not depend on the particular /-level under consideration, so that it is sufficient to introduce a single 
collisional relaxation rate for the whole term. As a consistency proof of our derivations, we showed that the transfer and relaxation 
rates for /-level populations and for interference between pairs of magnetic sublevels pertaining to the same /-level reduce to those 
derived in Sect. 7.13 of LL04 for the multilevel atom case. 

As an illustrative application, we considered a constant-property slab of given optical depth, located at a given height above 
the surface of a solar-like star, and anisotropically illuminated by its photospheric radiation field. The numerical solution of the 
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full non-LTE problem for the case of a two-term atom with unpolarized lower term shows that the polarization of the radiation 
emerging from the slab at ju = 0. 1 is sensitive to the presence of isotropic inelastic and superelastic collisions only for values of the 
parameter e' = CsifiuLuS — > P(L(S)/A(/3uLuS — > /3(L[S) on the order of 10"^ or larger. Such values are actually needed to produce 
an appreciable variation in the 2// profiles of the emergent radiation. 
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Appendix A: Properties of 3-7 and 6-7 symbols and of rotation matrices 

In this appendix we recall some useful properties and relations of 3- j and 6- j symbols, as well as of rotation matrices that are used 
in the derivation of the expressions presented in the paper A proof of these properties can be found in Chapter 2 of LL04. 

- Symmetry properties of 3- j symbols: 

The 3- j symbols are invariant under cyclic permutations of their columns and are multipUed by (-1)°+''+'^ under noncyclic ones 

a b c \ I b c a \ ^a+b+c I c b a 



R 1-1 « i = « r etc. (A.l) 

The 3-y' symbols are multiplied by (-1)°+*+'^ under sign inversion of the second row 

a b c\^^_^^,,,,J a b c \ ^^^2) 



\ a fi y j ' \ -a -/3 -j 

- Orthogonality relation of 3- j symbols 



»/3 

- Analytical expression of 3- j symbols for particular values of the arguments: 



= (-1)""" 5ah K-P ■ (A.4) 



- Symmetry properties of 6- j symbols: 

The 6- j symbols are invariant under interchange of any two columns and under interchange of the upper and lower arguments 
in any two columns. 

- Sum rules of 6- j symbols: 

+b+c+d+e+f+g+h+i+jnn ^ ^\^ " ^ cljfl b c \ i g h c \ _i f i j \ j f i j 



- Analytical expression of 6- j symbols for particular values of the arguments 

i ^ UrS..^.-,r-i r+^+/ ^ (A.7) 

\d e f j V(2fl+l)(2t/+l) 

- Contraction of 3-/' symbols: 

- Orthogonality relations of rotation matrices 

^ Di^iRyDi^iR) = 6mn . (A.9) 
p 

- Product of two rotation matrices: 

^)m^W^)m'^'W* = (-1)"Z^2/:+1)( ^ J^, ^J(jJ, ^, (A.io) 

%n(R)K'N'(R)-Y.(2K+ 1)1^1^ ^,^D'^Q,(Ry. (A.11) 
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